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Abstract 

We investigate the properties of the radion in the 5D hnear dilaton model arising 
ly-^ \ from Little String Theory. A Goldberger-Wise type mechanism is used to stabilise 

^f^ ' a large interbrane distance, with the dilaton now playing the role of the stabilising 

^^ ■ field. We consider the coupled fluctuations of the metric and dilaton fields and identify 

^^ ' the physical scalar modes of the system. The wavefunctions and masses of the radion 

and Kaluza-Klein modes are calculated, giving a radion mass of order the curvature 
scale. As a result of the direct coupling between the dilaton and Standard Model 
/\ , fields, the radion couples to the SM Lagrangian, in addition to the trace of the energy- 

momentum tensor. The effect of these additional interaction terms on the radion decay 
modes is investigated, with a notable increase in the branching fraction to 77. We 
also consider the effects of a non-minimal Higgs coupling to gravity, which introduces a 
mixing between the Higgs and radion modes. Finally, we calculate the production cross 
section of the radion at the LHC and use the current Higgs searches to place constraints 
on the parameter space. 
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1 Introduction 

For more than a decade extra dimensions have provided new insights into the hierarchy 
problem of the Standard Model (SM). By allowing fields to propagate in the extra dimensions 
the fundamental cutoff scales can be lowered. This was first considered in the Arkani-Hamed, 
Dimopoulos, Dvali (ADD) model p^, where if SM fields are confined to the boundary and 
only gravity propagates in the extra dimension, the fundamental scale of quantum gravity 
can be near the TeV scale. Assuming that the extra dimensions are stabilised, the weakness 
of gravity can then be attributed to the large volume occupied by the extra dimensions. An 
alternative possibility is to consider a warped extra dimension as in the Randall-Sundrum 
(RS) model \2\. A warp factor then controls the scales of the four-dimensional (4D) Minkowski 
slices, causing the local cutoff to depend on the location in the extra dimension. In this case 
the hierarchy problem is addressed by choosing the Higgs to have a TeV cutoff scale. This 
solution crucially relies on stabilising the extra dimension between the branes without any 
tuning of parameters. The Goldberger- Wise mechanism [3J , which introduces a bulk scalar 
field with boundary potentials, provides such a way to stabilise the interbrane distance. 

Recently a new solution generated by stacks of NS-5 branes in string theory was studied to 
also address the hierarchy problem. In the decoupling limit the NS-5 branes give rise to Little 
String Theory - a strongly-coupled nonlocal theory in six dimensions with no (apparent) 
Lagrangian description jHE]. However by holography [6j the seven-dimensional gravity dual 
description is a linear dilaton background which can be used to study the phenomenology of 
TeV Little String Theory [7j. An effective five-dimensional (5D) description can be obtained 
by compactifying two extra dimensions. This leads to the 5D linear dilaton model [8j which 
is compactified on a Z2 orbifold, with a 3-brane placed at each of the endpoints. The 5D 
linear dilaton model contains features of both the ADD and RS models, including a large 
volume bulk like in the ADD model, and a warp factor like in the RS model. Moreover, 
phenomenologically the graviton Kaluza-Klein (KK) spectrum is distinct, consisting of a 
~ TeV mass gap followed by a near continuum of KK resonances separated by only ~ 30 
GeV [8, 9J. 

However, just like in previous solutions to the hierarchy problem involving extra dimen- 
sions, the issue of stabilising the extra dimension is important. In this paper we investigate 
the stabilisation of the 5D linear dilaton model and study the resulting radion phenomenol- 
ogy. An interesting feature of this model is that it already comes equipped with a scalar field, 
the dilaton, that can be used to stabilise the extra dimension. By adding boundary poten- 
tials, as in the Goldberger- Wise mechanism, a stabilising potential can then be obtained. The 
usual simplifying assumption, following from Ref. [3], is to consider infinite boundary mass 
terms. However in our analysis we will relax this condition and wherever possible consider 
finite, or else finite but large boundary mass terms, similar to the general analysis considered 
in Ref [TOt [TT] . Apart from the usual radion couplings to the trace of the energy- momentum 
tensor, this leads to new radion couplings to SM fields that are confined to the brane. We 
present the Feynman rules and use them to study the radion phenomenology in this setup. 
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Similar to the graviton KK spectrum, the radion mass spectrum consists of a near con- 
tinuum of resonances spaced by ~ 30 GeV together with a massive radion zero mode. The 
radion zero-mode and lowest KK mode masses can be parametrically lighter than the ~ TeV 
mass gap, by suitably choosing the boundary mass parameters. The radion couples to the 
trace of the energy-momentum tensor in the usual way, but in addition there is a new direct 
coupling between the dilaton and the SM fields. The typical strength of these couplings is 
of order (10 TeV)~^, with the dilaton coupling being further suppressed by a factor inversely 
proportional to the boundary mass term. The SM fields have larger couplings to the radion 
compared to the near-continuum, and therefore the zero mode will be the first observable 
mode. 

Furthermore, as is well known, the radion can kinetically mix with the Higgs boson via 
a Higgs- curvature interaction, but since the typical coupling strength is of order (10 TeV)~^ 
there are negligible effects on Higgs phenomenology. However, the production and decay of 
the radion leads to observable signals at the LHC In particular due to the additional direct 
coupling between the radion and gauge boson kinetic terms, the branching fraction to 77 can 
be significantly enhanced. Since the decays are similar to that of the Higgs boson, current 
Higgs searches are then used to constrain the parameter space of the model. 

The outline of our paper is as follows. In Section 2 we briefly review the 5D linear dilaton 
model of Ref. [8]. In Section 3 we identify the radion and solve for the KK mass spectrum 
in the limit of large but finite boundary mass terms. The radion couplings to SM fields are 
computed in Section 4 where the effects of mixing from a Higgs-curvature interaction are also 
included. The phenomenology of the radion is presented in Section 5 where the radion decay 
widths, branching fractions and production at the LHC are studied. In particular recent 
LHC results are used to constrain the parameter space of the model. Finally in Section 6 we 
present our conclusion. Further details of the computations are included in the Appendices 
as well as the Feynman rules for the radion couplings. 

2 Linear Dilaton Model 

The set-up of the linear dilaton model is similar to the RS scenario, with a finite extra 
dimension, z, compactified on a Z2 orbifold, except that we will take the TeV scale to be 
the fundamental scale and the 4D Planck mass as the derived scale. There are two 3-branes 
located at z = and z = r^ with the SM confined to the z = brane. We also introduce 
a 5D bulk scalar field known as the dilaton. In the Einstein frame we have the following 
action [8] 

Sbulk = / d^xy/^ 



M''(Tl--{y'^f\ -ei^A 



(2.1) 



Svis{hid) — d x\/ —ge'A'^ \l^vis{hid) — K)is(hjd)) , (2.2) 



where i^ is the dilaton field and M is the fundamental scale, which is of order the TeV scale. 
In the conformal coordinate, z, the background solutions to the field equations are given by 

(^{z) = a\z\, (2.3) 

ds^ = e-i°l^l {ri^^dx^'dx" + dz^) , (2.4) 

where (f){z) denotes the background value of the dilaton field, </?, and |a| < ^4= is required to 
ensure the 5D curvature is smaller than the fundamental scale. There are also the following 
constraints 

A = -M^a\ V^is = -Vhid = ^aM\ (2.5) 

which correspond to the usual tuning of the 4D cosmological constant. 

Taking the curvature term in the action (12.11) and integrating over the fifth dimension 
yields the effective 4D Planck mass, 

dz e-°l"lM3 = -2 (e"""^ - l) . (2.6) 

We then clearly require that a < and lard ~ 70 in order to obtain the required value for 
Mpi. 

In order to allow a straightforward comparison with the RS model, we can perform the 
following coordinate transformation dy = e's'^^dz. In the new coordinate, the metric be- 
comes 

ds"^ = (l + \^ ) r^^^dx^'dx" + dy"^. (2.7) 



We can now see that this model exhibits a power-law warping as opposed to the exponential 
warping of the RS model. The size of the extra dimension in this case must therefore be 
substantially larger than in the RS model in order to obtain the necessary hierarchy. Taking 
|a| ~ TeV and enforcing the correct value for Mpi allows us to determine the proper length 
of the extra dimension, giving y^ ~ 10 nm. This compares with a proper length of ~ 0.1 mm 
in the ADD model (for two extra dimensions) and ~ 10~^^ cm in the RSI model. 

The linear dilaton model has an interesting spectrum of graviton KK modes. There is a 
single massless mode, which has a flat profile in the extra dimension and is identified with the 
usual 4D graviton. In addition, there is a KK tower of excited states, with a mass spectrum 
given by 

^n = -r+[ — ] , (2.8) 
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where n = 1, 2, 3, ... . Of particular interest is a large mass gap between the massless graviton 
and the first of the KK modes, which are then closely spaced. The mass gap is determined 
by the curvature scale, a, and for a TeV mass gap we find that Vc ~ (30 GeV)^^ and the KK 
modes essentially form a continuum of states. The graviton KK modes are localised near the 



SM {z = 0) brane, and their couplings to the Standard Model fields are given in Ref. [8]. 
They are found to couple to the Standard Model fields with a strength of order (80 TeV)~^ 
for the lowest mode. 

2.1 Stabilisation 

In order to truly solve the hierarchy problem the interbrane distance must be stabilised. In 
the linear dilaton solution of [8J this can be done by using a variation of the Goldberger- 
Wise mechanism. Unlike the Randall-Sundrum scenario, the linear dilaton model naturally 
includes a bulk scalar field, the dilaton, that can in fact play the role of the stabilising field 
once it acquires a vacuum expectation value on the branes. 
We begin by taking boundary potentials of the form 

Vvis(hid){^) = ^Kis(hid) + l^vis{hid)M {if — (f)vis(hid)) , (2.9) 

where Xvis{hid), fJ'vis{hid) and (f)vis{hid) are constants. At this point it is useful to formulate the 
model in terms of the solution generating method of [12], by considering the superpotential 
W{(f) = WqC^ , where Wq is a constant. The background field equations include two junction 
conditions at the boundaries, which can now be expressed as 

iy(0) = ±elv,i,ih^dM, (2.10) 

Using (12. 9p . these equations can only be simultaneously satisfied if = (pmsihid) on the 
boundaries. Thus we see that (p has developed a non-zero vacuum expectation value on 
the branes as a result of imposing the boundary potentials (12. 9p (setting fims{hid) = also 
satisfies the conditions, however this simply reproduces the unstabilised case). It is worth 
noting that an additional linear term can be added to the boundary potential, however this 
simply shifts the value of the VEV and the potential can always be rewritten in the form 
(12. 9p . The requirement that = (f)ms{hid) on the boundaries then gives 

Xvisihid) = Wo = aV^AAP, (2.12) 

where we have used the definition of the bulk potential from the action (12. ip . The condition 
(I2.12P corresponds to the usual tuning of the 4D cosmological constant in (12. 5p necessary to 
obtain a fiat brane solution. 

The background solution for the bulk scalar is once again given by 

(j){z) = a\z\ + (f), (2.13) 

where here we have included the integration constant, 0, which was previously taken to be 
zero. Imposing the boundary conditions = (p^s and = (phid on the z = and z = re 



branes respectively, we then find that 

= (pvis, (2.14) 

arc = (phid - (pvis- (2.15) 

The second expression clearly shows that we have indeed stabilised the brane separation, r^. 
In order to achieve the correct hierarchy, we previously required that lard ~ '''O- This can 
be achieved without the need for any extreme fine-tuning of the parameters. In the proper 
coordinates ( 12. 7 P this corresponds to stabilising the interbrane distance at y^ ~ 10^°|a|^^, 
much larger than in RS models. Finally, it is important to note that since ipms appears in the 
exponential dilaton factor in the boundary action (12. 2p . we require that 0^js = in order to 
successfully reproduce the Standard Model. This is the case we shall consider from now on. 

3 Identification of the Radion 

In this section we identify the physical radion and KK modes and calculate the mass spectrum. 
The scalar fluctuations about the background solutions are denoted 

ds^ = e-i"l"l [(1 + 2^')ri^^dx^'dx'' + (1 + 2^)dz^] , (3.1) 

(^(x,z) = 0(z) + (50(x,z). (3.2) 

We then obtain the linearised Einstein equations for the perturbations, which include the 
following constraint equations 

$ + -vl/' + hep = 0, 

a 3 (3.3) 

$ = -2^, 

where primes denote differentiation with respect to z. The presence of these non-dynamical 
equations reflects the fact that the gravitational and dilaton scalar perturbations are coupled 
and the action must be diagonalised in order to identify the physical degree of freedom of 
the system. Using these constraint equations the dynamical equation can be expressed as 



6-5"^$) =0. (3.4) 



d^ 0?^ 

There are also two junction conditions at the boundaries. The first is obtained by simply 
evaluating (13. 3 P at the boundaries, while the second contains additional physical information 
and is given by 

4M^e^"^ 

H = ±-^7^^ ^ (^'^ - "*) ' (3.5) 

g^ I 3e 3 Vyis(hid) (0) 



One of the nice features of this model is that the equations of motion for the perturbations 
can be solved analytically. While the boundary conditions are non-trivial, they become 
considerably simpler in the limit -^{e^V) — )■ oo, which forces 50 — )■ on the branes. 
This is the limit most often considered in the literature, including the original Goldberger- 
Wise case. However, this will be insufficient for our purposes since many of the interesting 
phenomenological features of this model arise due to the non-zero coupling of the dilaton 
field to the SM Lagrangian on the brane, as shall be seen in sections H] and IS 
We begin by performing a KK decomposition 



$(x^z) = 5;]$„(^)Q, 



\x'^ 



(3.6) 



n=0 



where Qn{x^) are the 4D Kaluza-Klein modes and satisfy the Klein-Gordon equation, and 
$„ are the profiles in the fifth dimension. The equation of motion is then simply given by 



dz^ 4 



with the boundary condition 



2a " ' 






A$" 



2a ± 9/i„is(/jirf) \2a 
We obtain the following solution for $„, 



0, 



3$1 - «$, 



^r 



N„e-^' 



sm{(3nZ 



■ COs{l3nZ) 



(3.7) 



(3.8) 



(3.9) 



where we have defined /3^ 



mz 



a 
4 ' 



4/32 + a{a- fx^is 
Nr, is an overall normalisation factor, and we have 



used the boundary condition at z = to fix one of the constants. The mass spectrum is then 
obtained by evaluating the boundary condition at the z = r^ boundary. The complete mass 
spectrum can be determined analytically by taking ^vis{hid) to be large and solving as an 
expansion in \a\/ ^ms(hid)- However, an exact expression can also be obtained for the radion 
and lowest KK mode in the limit \arc\ ^ 1. We then find that 

2 



m. 



rad 



a 

T 



a 



IQei 



3-A/9 + 4e^i, + 4e 



< e^is < cx). 



(3.10) 



where eyis(hid) = |«|//^«s(/ijd) and fivis{hid) > is required to ensure there are no tachyonic 
modes. Note that dependence on thid in (I3.10p is suppressed in the {ard ^ 1 limit. We also 
solve the boundary condition as an expansion in eyis(hid) to obtain the complete KK spectrum 
at first order 



"^rad 



m^ 



-a 
9 



1- 



ra^TT^ 



2e, 



4(4nV + 3|aren 
12'n?T['^\arc\ + \aTcV' 



+ ^hid) 



(3.11) 
(3.12) 



where n = 1,2,3, ... . Note that in the hmit ems{hid) -^ 0, corresponding to fims{hid) -^ oo, we 
obtain the resuh given in [lOj. 

We see that this model also gives rise to a rather interesting mass spectrum for the scalar 
fluctuations. Both the radion and lowest KK mode have masses of order the curvature scale, 
but can be parametrically lighter depending on the values of ems{hid) ■ The spacing of the KK 
modes is then determined by the size of the extra dimension as we would expect. Once again 
we have Tc ~ (30 GeV)~^ and the KK spectrum forms a near continuum of modes above a 
mass of approximately half the curvature scale. 

The mass spectrum is plotted as a function of e in Figure [H where the KK modes have 
been determined numerically and we have taken e^g = etid- Note however that as seen in 
(I3.10p the lowest mode is dependent only on eyis, while the second lowest mode depends only 
on thid for Chid ^ 1- We find that the exact and approximate solutions are well matched 
when e is small, justifying our expansion. Interestingly, we observe that both the radion 
and lowest KK mode are highly sensitive to the boundary mass terms and as e is increased 
there are two distinct modes below the near continuum of KK states. In the e — ?■ oo limit, 
corresponding to the unstabilised case with fims(hid) = 0, the KK spectrum is identical to the 
graviton spectrum f|2.8p with a ~ TeV mass gap. We note that while the decomposition (13.1 p 
becomes ambiguous when considering massless scalar modes, we expect two massless modes 
in this limit, where the additional massless mode arises due to the extra bulk scalar field (the 
dilaton) . 
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Figure 1: The radion and KK mass spectrum as a function of e, where we have taken 
^vis = ^hid- The lowest five modes are shown and the dashed lines denote the approximate 
solutions. 



4 Coupling to the Standard Model 

In this section we calculate the couplings of the radion and KK modes to the Standard Model 
fields. The relevant term in the boundary action is 



S = jd^x^-ge-CsM. (4.1) 

where Csm is the SM Lagrangian and g is the induced metric on the z = brane. We clearly 
see that the radion couples to the SM fields through the induced metric on the brane and also 
via the exponential dilaton factor. This dilaton factor is not present in the usual RS models 
and gives rise to an additional coupling between the radion and SM fields. Expanding the 
above action to first order in the perturbations, we obtain the following interaction action 
between the radion KK modes, Q^, and the SM fields 

S^nt=\Y,^n{'^) j d^X^/^0Qn9^''T^u + \Y,5(l)n{'^) I d^X^/^0QnCsM\^^^^^ (4.2) 

n n 

where ^o is the induced background metric on the brane and T^j, is the Standard Model 
energy-momentum tensor. The first term in (14.21) is the usual coupling of the radion to the 
trace of the energy momentum tensor, while the second term arises due to the presence of 
the dilaton and produces a direct coupling between the radion and SM Lagrangian. This 
additional coupling can have a significant effect on the radion phenomenology, as we shall 
see in section |5l 

For convenience we define the following coupling constants 

l^^,n ^ '^n(O) «^ ^ ^0n(O) ,. „s 

M ~ 2 ' M ~ 3 ■ ^ ' ' 

The couplings can now be determined by correctly normalising the solutions for the pertur- 
bations obtained in section [31 As we found in section [HI the couplings for the radion can be 
determined analytically in the jard ^ 1 limit, while the KK couplings can only be obtained 
as an expansion in ems(hid)- The exact expression for the radion coupling is non-trivial and 
does not provide additional insight, so here we give only the approximate solutions. The 
couplings to the energy-momentum tensor are given by 

\l^<i>rad\ = ^A/^TT IH ^ots 5 (4.4) 




^^,n\ = ^^, r;77\/ T7 (1 ~ ^«^) ' (4-5) 



where we have also taken \aT^ ^ 1, since the expressions simplify significantly in this limit. 
We see that the coupling strength is largely determined by the fundamental scale, M and 



also the ratio a/M. Additionally, the couplings for the KK modes are suppressed by a factor 
of I arc I ^/^ relative to the radion mode. Hence, we always expect the single radion mode to 
be observed first by experiments, even when the radion mass is close to that of the KK near 
continuum or when there are two light modes. It is also interesting to note that the couplings 
for the KK modes increase with n, however this remains true only for small n since fl4.5p is 
no longer valid when n^ ~ la'^d- 

brane 

(4.6) 
(4.7) 



Using the constraint equation for the perturbations (13.31) evaluated at the z 
we can also determine the couplings to the SM Lagrangian, which are given by 



\l^(j),rad\ 
\^it>,n\ 




We see that unlike nq, above, these couplings do not include a zeroth order term, consistent 
with the fact that 50 — )• in the /i^j^ — )• oo limit. The couplings to the SM Lagrangian 
are therefore suppressed by an additional factor of \oi\/ ix^s) which must be small in order 
to justify our expansion. Once again, the KK modes are further suppressed relative to the 
radion mode. 
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Figure 2: The radion couplings as a function of e^is. The upper and lower curves correspond 
to «;$ and k,^ respectively. The dashed lines denote the approximate solutions. 



Figure [2] shows the radion couplings as a function of ems-, where the upper and lower 
curves correspond to /«$ and k^ respectively. Once again we see that the approximate and 
exact solutions are in good agreement when e^s is small. We also find that /«$ > k^. 
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even for larger values of e^s where our approximate solution is no longer valid. Taking 
M ~ TeV, |a|/M ~ 0.5 and ems = 3, we find that K<s,^rad ~ 0.11 and K,<s,^rad/M will be of order 
~ (9 TeV)""*^. Similarly, K^^rad ~ 0.07 and K^^rad/M will be of order ~ (14 TeV)~"^. Note also 
that while the overall sign of the couplings is undetermined by the normalisation condition, 
the relative sign is fixed giving k^/k,^ > 0. 

4.1 Higgs-curvature interaction 

An interesting situation arises for scalar fields propagating on the brane, since they can 
couple non-minimally to gravity. This allows us to introduce the following additional Higgs- 
curvature interaction term to the 4D effective action. 

S^= [(fx^^e^in{g)H^H, (4.8) 

where 'R-{g) is the Ricci scalar obtained from the 4D induced metric on the brane. For 
an arbitrary scalar field if), we could also in principle include the linear coupling 7l{g)ijj, 
however in the case of the Higgs this is forbidden by the gauge symmetry. As we shall see, 
this additional interaction term results in a mixing between the Higgs and the radion. The 
analysis given below follows that performed for the RS case in [T3[ IT4] . however here we 
include the additional effects of the dilaton field. 

The 4D Ricci scalar is calculated from the induced metric to be 

n{g) = -6(1 - $)-3/V"5A(l - '^)'/', (4.9) 

which results in the following Lagrangian for the interaction 

C^ = -6^6^(1 - <^Y/^ri>"'d^d,il - ^f'^H^H. (4.10) 

Now, after expanding the Higgs field around its vacuum expectation value, H = {h + v)/\/2, 
where v ~ 246 GeV, to quadratic order in the fields we have 

C^ = -3^ (jv^<^D<^ - -v^6cf>D(l> - vhD<^] . (4.11) 

Finally, inserting the KK expansion for the fields and considering only the radion mode, r, 
we obtain the following complete Higgs-radion Lagrangian 



C = — hUh- -mlh'^ - - 
2 2 '^ 2 



M2 ^'^^ ~ "^^^ 



rDr - -rn^y + -^^hDr, (4.12) 



where ruh, rrir are the Higgs and radion masses respectively, in the ^ = limit. The ^ term 
then clearly introduces a kinetic mixing between the Higgs and radion fields. The Lagrangian 
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must therefore be diagonalised to identify the physical mass eigenstates of the system. We 
begin by diagonahsing the kinetic terms via the following transformation 



h = h -\ — r 



where 



n' 



1 + 



M2 



O' 



1 -6^)k^ -«</,). 



(4.13) 



(4.14) 



This in fact also enables us to place a constraint on the parameter ^, since Q"^ must be 
positive in order to ensure that the radion mass term retains the correct sign and we avoid 
encountering a ghost radion field. We therefore require that C, lies in the range 



12k<s, 



P 



o 4M2 \ ^ 1 



12/t$ 



P+\ P' + 



4M2 



(4.15) 



where p = k,,^ — k^. Taking k$ = 0.11, k^ = 0.07 and M = 1 TeV, we find that this places a 
limit of-6.1 < ^< 6.2. 

So far, we have successfully diagonalised the kinetic terms in the Lagrangian, however 
this now introduces a mixing in the mass matrix, which must be diagonalised via the rotation 



h' = cos Ohm + sin^r^, 
r' = — sin 9hm + cos 9rm, 

where hm, Vm are the physical mass eigenstates and the mixing angle is given by 



(4.16) 



tan 261 



12^K,q>vQ 



mt 



M 



ml — m\ ( 1)2 






(4.17) 



We can now express the gauge eigenstates in terms of the mass eigenstates as 
h-- 



cos U sm o \ nm+ \ sm -\ — ^ ^ ^ cos U r^. 



QM 



QM 



sin 6 



hr 



COS I 



(4.18) 



We can also calculate the corresponding mass eigenvalues, which are given by 

1 



m^ 



29? 



m, + m^ I 1 + 



M2 



± 



m^ — m 



6C/t|t;2 
M2 



2 1 ^ ^'-2,^2„,2^2^2^ ' 



+ 



144C, K^v^m^mf^ 
M2 



(4.19) 
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In the following section we shall only be considering the radion interactions described in 
(14. 2 p and will not be concerned with higher- dimensional operators, which are suppressed by 
higher powers of M. It therefore makes sense to expand the above results to leading order 
in 1/M. The mixing angle then becomes 

tan 2^ ^ —^ ^ ^ ^ . 4.20 

Now, provided that mr^rrih, the mixing angle is small and our results simplify to 



'^ (4.21) 



M mi — mi 



QiK^v m^ 

''"m ' I m- 



M ml — mj\ 



This mixing can in principle have an interesting effect on the Higgs phenomenology, since 
a strong mixing would lead to a reduced cross section for the Higgs. However, ( 14.2ip shows 
that the extent of the mixing is largely determined by the factor hiq,v/M and is therefore 
small, except in the case when mr is very near to mh-, or if we take ^ to be large. While the 
mixing may not have any significant effect on the Higgs, it still plays an important role in 
determining the radion phenomenology, since the extent of the mixing is comparable to the 
size of the radion couplings. 

When m^r^m^h, (I4.2ip is no longer valid and we must include higher order terms in the 
expression for the mixing angle (I4.20p or else consider the full expressions in (I4.18p . In this 
case we find that the mixing becomes large and the two mass eigenstates are now essentially 
part Higgs and part radion. The result is two states with closely spaced masses, both of 
which interact in a similar way to the SM Higgs but with a reduced cross section. 

We also note that at leading order our result does not depend on k,^ and the dilaton does 
not contribute significantly to the Higgs-radion mixing. Our results are then equivalent to 
those for the RS case given in [T3| [l4j. In this limit, the mass eigenvalues also simply reduce 
to m,r, m,h. The above analysis can also be performed with the inclusion of the radion KK 
modes and is given in Appendix [Bl The mixing effect is sub-dominant compared to the radion 
only case. 

4.2 SM-Radion interactions 

As we have shown, the radion couples to the trace of the Standard Model energy-momentum 
tensor as well as directly to the SM Lagrangian. We now derive in detail the radion interac- 
tions with the Standard Model fields and present the corresponding Feynman rules for the 
interaction terms in Appendix ICl 

Following on from our discussion of the previous section, we consider the case where there 
is a mixing between the radion and Higgs fields and express the interactions in terms of the 
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physical mass eigenstates. In order to simplify the expressions, we shall use the following 
general formulae relating the gauge and mass eigenstates, where the values of the coefficients 
can be read off from fl4.18p : 

h = aohm + airm, ,^ ^^, 

r = bohm + hrm. 

We shall also restrict ourselves to interactions up to quadratic order in the SM fields, since 
these will be the dominant processes when considering the radion decay modes in the following 
section. The Standard Model energy-momentum tensor is given by 



3 

2 



+ ^ -9^(^*7'^$) - 3^i7'^9^^ + Am^^^. (4.23) 
/ 

We begin by considering the interactions of the radion with the massive gauge bosons, 
which are given in the unitary gauge by 

Ant = -^r^ {2WlW^'^ + ^..^^'1 + (^(^ - ^'^) + ^) rm i^ml^WlW^ + m\Z,Z^) . 

(4.24) 
We see that as a result of the dilaton coupling to the SM Lagrangian, there is now a coupling 
between the radion and the gauge boson kinetic terms, which will become important at large 
momenta. As we shall see, this additional coupling is particularly significant for the massless 
gauge fields. 

We now consider the interactions with the SM fermions which are given by 

(4.25) 
However, in most of the cases we shall be considering it will be sufficient to restrict ourselves 
to on-shell fermions. Then using the Dirac equation we can significantly simplify the above 
interaction and obtain 

Ant = m^ (^ - ^) r^^M/. (4.26) 

Of course, the Dirac Lagrangian vanishes after using the equations of motion and our inter- 
action term is therefore independent of k^. 

There are additional considerations when looking at the radion interactions with the 
Higgs boson. As a result of the mixing between the radion and the Higgs, the interactions 
will now be dependent on the radion potential as well as the Higgs self-interaction terms. In 
addition to this, there are interaction terms which arise from the Higgs-curvature term (14.81) . 
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Following the analysis in section 14.11 but here considering terms which are cubic in the fields, 
gives the following interaction term to leading order in 1/M 

>CL = ^^^^^ {d.kmd^h^ + h^Uh^) r„. (4.27) 

Now putting this together with the other interaction terms we obtain 



'~-int — fi/T \ ^^ ijK<J> + yrn(J^j,i^mC) l>,m \ , , Tm^-mLJ/l^ 



2 2 / 1 /n ^<t>\ iOil 



+ aX^^(2«:*-f)-^J^r./^™, (4.28) 

where the final term arises from the cubic Higgs self-interaction term in the SM Lagrangian. 
We have excluded the effects of the radion potential terms, which we shall assume to be 
small. 

Finally we consider the interactions with the massless gauge bosons. Here we must also 
include the quantum effects of the scale anomaly, which is reviewed in [15]. The interaction 
with the gluons is given by 

£.„. = "5 (^ + f!Ag£^^ ,^^._^^„._ (4.23) 

where hqcD = 11 — 2ri//3 with nj = 6. For the gluons it turns out that the direct coupling of 
the dilaton and the additional coupling due to the scale anomaly are of comparable strength 
due to the relatively large value of a^ and the suppression of k,^ relative to k$. 
The interaction with the photon is similarly 



An. = -ii 7 + ^^T^ih + by) r^F,^F>^^, (4.30) 



bi /K^ OtEM^^'^ 

where 62 = 19/6 and by = —41/6. However in this case the anomalous contributions to T/^ 
are insignificant due to the much smaller value of aEM- We similarly ignored the effects of 
the scale anomaly when considering the W and Z gauge bosons above. 

5 Radion Decays and Production 

5.1 Decay widths 

From the interaction terms given above and the Feynman rules given in Appendix [C] we 
calculate the partial decay widths of the radion into WW, ZZ, hh, ff, 77, and gg. We have 
taken our results to leading order in 1/M and used the expression 04.2ip for the Higgs-radion 
mixing. They are given by 
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T(r^W^W~ 



ra" 



167rM2 



a/1 — Xw 



-^ I" /i^^ I 1 ^ Xw + T^PF 



6^m^ 



1- 
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mt — m 



-K(f)K,^ ( 1 + —Xw 
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1, 



-rfr^py+w^- 



(5.1) 



(5.2) 
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327rM2 
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Xh 
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(5.3) 



(5.4) 
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87rM2 
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27r 
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m^ — m 



(5.5) 



where 



/(:^) 



sm 

1 

' 4 



n 2 



, 1 + yr^^ . 

In , — m 



1 2 



x> 1, 



X <1, 



(5.6) 



(5.7) 



and we have defined Xi = Amj/m,"^, {i = W,Z,h,f) and Nc = 3(1) for quarks (leptons). 
When calculating the partial decay widths to gg and 77, we have also included the one-loop 
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contributions of the top quark and the W boson, which become important when the mixing 
is strong. We have neglected one-loop corrections proportional to k^. Finally, in the limit 
i = and K0 — )■ our results reduce to those for the RS case given in [161 HT] - 

5.2 Branching fractions 

In this section we use our previous results to look at the branching fractions of the various 
radion decay modes. We begin by considering the case where there is no mixing between the 
radion and the Higgs (^ = 0). Figure [3] shows the radion branching fractions as a function of 
its mass where we have taken e^g = 3 and M = 4 TeV. We have chosen a value for the Higgs 
mass of 125 GeV. 

We see that in the low mass range the branching fractions are dominated by the gg and 
77 channels, while in the high mass range gg and WW are the dominant decay modes, but 
with significant contributions from several other channels. The large branching fraction to 
77 is of particular interest and is a direct result of the dilaton coupling to the gauge boson 
kinetic terms. It also remains significant even at high radion masses. 

If we compare the branching fractions for the linear dilaton model shown here with those 
of the RS model [16j, we note that the most important difference comes in the branching 
fraction to 77. This is significantly enhanced in the linear dilaton model, since the dilaton 
coupling is much stronger than the scale anomaly and one-loop contributions, which are 
proportional to a em- We conclude that the large branching fraction to photons seen in 
figure E] is therefore an important distinguishing feature of this model. 

Next we look at the case where there is a mixing between the radion and Higgs. We 
shall take C, = 1/6, which corresponds to the conformal limit when we also take m/i=0. The 
radion branching fractions are then shown in figure HJ We observe that at high mass, the 
situation is similar to the zero mixing case in figure [21 except that the hh decay mode is now 
suppressed as a result of the mixing. In the low mass region we observe an interesting drop 
in the branching fractions to gg and 77 when the radion mass is near to that of the Higgs. 
This is the result of a cancellation in the partial widths (15. 5p . (15. 6p due to the strong mixing 
in this region. 
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Figure 3: Branching fractions of r^ as a function of its mass with e^i^ = 3, M = 4 TeV, ■C = 
and ruh = 125 GeV. The left and right panels are the same but cover a different range in 
mass. 
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Figure 4: Branching fractions of r^ as a function of its mass with e^g = 3, M = 4 TeV, 
^ = 1/6 and ruh = 125 GeV. The left and right panels are the same but cover a different 
range in mass. 



5.3 Production 

In this section we briefly discuss the production of the radion at colliders. At hadronic 
colliders, such as the LHC, radion production will be dominated by gluon-gluon fusion. The 
production cross section for this process at a hadron collider with centre of mass energy ^/s 
is given by 



a{pp-^r) 



dx , , (ml 



-g{x,mr)g 



\ sx 



TTlr 



mt 



^{gg^r) 



{5.i 



where g{x,q) is the gluon parton distribution function at momentum fraction x and renor- 
malisation scale q, and a{gg^r) is the gluon-gluon fusion subprocess cross section, which is 
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given by 



^{gg^r) 



TX 



64M2 



H + 17-^^ ( ^QCD + a;* (1 + (1 - xt)f{xt)) j 



(5.9) 



The gluon-gluon fusion production cross section at the LHC for a centre of mass energy of 8 
TeV is shown in figure [5] for several values of the fundamental sale, M. The parameterisation 
of CTEQ5L [18j has been used for the parton distribution functions. We find that even for 
relatively large values of M the radion production cross section remains significant. This 
suggests good prospects for searches at the LHC, as we shall see in the following section. 
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Figure 5: Gluon-gluon fusion production cross section for the radion at the LHC with a 
centre of mass energy of 8 TeV. We have taken e^s = 3 and ^ = 0. 



5.4 Constraints on parameter space 

Similarities between the radion and the Higgs boson enable us to place constraints on the 
parameters of the model using the results from Higgs searches currently underway at the LHC. 
Similar studies have recently been performed for the RS radion in [19] [20| |2T| [22| [231 |2l]. We 
begin by using a Breit-Wigner narrow width approximation to relate the radion and Higgs 
cross sections according to 



T{r^gg) Br{r-^X) 
asMipp^h^X) ^ r{h^gg) Br{h^Xy 



(5.10) 



This ratio is plotted as a function of the radion mass for the WW, ZZ and 77 decay 
channels in figure [6] We observe that the strongest constraints on the model will be obtained 



using the 77 decay channel, which gives a value of a josu of order one, due to the large 
branching fraction of the radion to this mode. 
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Figure 6: o Josm as a function of the radion mass. We have taken e^s = 3, M = 4 TeV and 



We use the latest results from the ATLAS experiment in the /i — )■ 77 channel with 4.9 fb^^ 
integrated luminosity [25], where the CL^ method has been used to place bounds on o Josm 
at the 95% confidence level. Comparing these results to (I5.10p we can place constraints on 
the parameters of the model, namely the fundamental scale, M, and the curvature scale, 
a. These constraints are shown in figure [TJ where the lighter shaded region is excluded by 
the data, and the darker region is the theoretically excluded region where the 5D curvature, 
'^o? is larger than M^. The dashed line shows the constraints from ADD total cross section 
searches at the LHC [S]. 

We obtain a limit of M > 3 TeV for values of \a\ in the range 340-460 GeV. Note 
however that these limits will depend on the value taken for e„js, which is not constrained. 
For larger values of e^s the limits become stronger, while for e^s ^ 3 the ADD total cross 
section searches provide the best constraints on a and M. Using the WW ^Ivlv and ZZ^Al 
channels, constraints can also be placed on the parameter space in the higher mass regions. 
However, due to the smaller value of (t/ctsm for these channels, the limits currently lie below 
the theoretical exclusion. 
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Figure 7: Constraints on the model parameters, a and M, using ATLAS results in the 
/i — )■ 77 channel. The lighter shaded region is excluded by the data and the darker region is 



the theoretically excluded region where M^ < 
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The dashed line shows the constraints 



from ADD total cross section searches at the LHC |9j. We have taken e^s = 3. 



6 Conclusion 

Motivated by Little String Theory, the 5D linear dilaton model provides a distinctive solution 
to the hierarchy problem. It is characterised by a graviton KK mass spectrum which consists 
of a ~ TeV scale mass gap, followed by a near continuum of resonances spaced ~ 30 GeV 
apart. The scalar sector of the model plays a crucial role, since it contains a scalar field, 
the dilaton, which can be used to stabilise the extra dimension. Just like in the Goldberger- 
Wise mechanism, the dilaton acquires a non-zero vacuum expectation value on the branes 
after imposing appropriate boundary potentials. A large volume is then stabilised with an 



10 



10 



a\ 



The correct hierarchy is achieved without the 



interbrane proper distance of yc 
need for any extreme fine tuning. 

The coupled scalar fluctuations about the metric and dilaton flelds were then considered. 
We solved the equations of motion for the system and obtained the wavefunctions and masses 
of the modes. It was found that the radion and lowest KK mode both obtain masses of 
order the curvature scale, but can be parametrically lighter depending on the boundary mass 
parameters. The KK modes were then closely spaced and essentially formed a near continuum 
of modes of order the curvature scale and above. 

This model introduces a direct coupling in the boundary action between the dilaton and 
the Standard Model flelds. It was shown that this results in a coupling between the radion and 
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SM Lagrangian, in addition to the standard coupling to the trace of the energy-momentum 
tensor. The strength of these couphngs was calculated, with a typical strength of order 
(10 TeV)~^, and a suppression in the coupling to the SM Lagrangian by a factor inversely 
proportional to the boundary mass term. Additionally, the couplings of the KK modes were 
suppressed relative to the radion, so that we would initially expect to observe only a single 
mode in experiments, despite the closely spaced mass spectrum. 

We also considered the situation where the Higgs couples non-minimally to gravity and 
included an additional Higgs-curvature interaction term in the 4D effective action, resulting 
in a kinetic mixing between the radion and the Higgs. However, we found that at leading 
order the dilaton did not contribute to the mixing and our results reduced to those for the 
RS case. The mixing then had virtually no effect on the Higgs phenomenology, except in the 
case when the difference between the radion and Higgs masses was small. The contribution 
of the KK modes to the mixing was also investigated. However, we found that despite their 
close mass spacings, the cumulative effect of the KK modes was still too small to significantly 
effect the Higgs phenomenology. 

The decay widths and branching fractions of the radion were calculated in section \5\ and 
found to be substantially different to the RS radion. This was largely due to the additional 
coupling between the radion and the gauge boson kinetic terms. In particular this led to a 
significantly increased branching fraction to 77, which provides an interesting distinguishing 
feature of this model. 

The production of the radion via gluon-gluon fusion at the LHC was also discussed and 
found to be significant even for relatively large values of M. Finally, we demonstrated the 
potential of current Higgs searches at the LHC to constrain the parameter space of this model. 
Using the latest ATLAS results in the /i — > 77 channel we obtained a limit of M > 3 TeV 
for values of \a\ in the range 340-460 GeV when fi^s = \ 
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A Wavefunction Normalisation 

In order to calculate the couplings to the Standard Model fields we need to ensure that the Qn 
modes are canonically normalised, which corresponds to correctly setting the normalisation 
constant A^„. However, the equations of motion are only linear in the perturbations and do 
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not contain sufficient information to correctly normalise the modes. It therefore becomes 
necessary to diagonalise the action (12. ip to second order in the perturbations. This analysis 
has been done for the general case in |20j. Using their result, we obtain the following free 
action for the 4D physical modes of the system, 



St' = C^ j d'^xQnp- ml]Qr,, 



27M3 r- , 
Cn = ^^l dze 



n 4a ^_ , 2a ^^ 



" 3 " " ^ 3 



-$. 



1 
2' 



(A.l) 



(A.2) 



In addition to correctly normalising the modes, diagonalising the action to second order 
in the perturbations allows us to determine the precise linear combination of the fields $, 50, 
which corresponds to the 5D dynamical variable in the bulk. 



X 



^^-"(l-f)' 



(A.3) 



where x is the dynamical bulk field and obeys the following equation of motion 



- d^ a2- 

dz^ 4 



x = o. 



(A.4) 



Here we provide the complete normalised wavefunctions for the radion and KK modes. 
The wavefunctions take the form of (13. 9p . where the normalisation constant Nn is given by 



N„, 



3M3/2 



Al3lae + a^{l + e) 



6/3„r,(4/3^ + a'){16/3y + 4/3>^(9 + 2e + 2e') + a\l + ey) 



- 8/3„a ( 16/3^e(9 + e) + 8/3>^e(l + e) + a^(-8 - 7e + e^) J sin(/3„rj^ 
+ 3(64/3^ + 16/3^^-9 - 2e + e") - 4/3^a'(-6 + 4e + e") - a^(l + e)A sin(2/3„re) 



A/2 



(A.5) 



where we have defined e = e^g for simplicity and /3„ is defined near (13. 9p and can depend on 
both e^is and ehid- 
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B Higgs-radion KK mixing 



We consider the effect of the Higgs mixing with the closely spaced radion KK modes. The 
Lagrangian that needs to be diagonalised is now given by 

c = -hah - \mlh^ + -£ (-iQ„nQ„ - \mlQl + ^^/^ng.) 



m n 

We see that in addition to the Higgs-KK mixing, the last term also introduces a kinetic 
mixing between the KK modes. However, this term is suppressed by an additional factor 
of 1/M. We shall therefore work to leading order in 1/M, which allows us to neglect the 
mixing between the KK modes and focus instead on the mixing with the Higgs. The analysis 
can now proceed as for the single radion case and we diagonalise the kinetic terms via the 
transformation 

h = h' + ^Yl ^-^Mn- (B.2) 

n 

The Lagrangian now takes the form 

£ = -\uuh! - i 5^g„nQ„ - i/wv^ (B.3) 



where we have defined the vector 



and the mass matrix is given by 



m*-' = vfii. 
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(B.4) 



(B.5) 



It is now straightforward to diagonalise the mass matrix numerically. We find that the 
mixing with the Higgs remains small and there is one mass eigenstate which is largely com- 
posed of the Higgs gauge eigenstate. So despite their relatively close mass spacing, the 
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cumulative effect of the KK modes is still too small to have a significant effect on the Higgs 
phenomenology. There is however an exception when one of the KK modes has a mass very 
close to the mass of the Higgs, in which case the mixing between those two states becomes 
significant. However the mixing with the other KK modes reduces this effect compared to 
the radion only case discussed previously. 

C Feynman Rules 
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